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ABSTRACT 

The metallic gas associated with the /3 Pic debris disk is not believed to be primordial, but arises from 
the destruction of dust grains. Recent observations have shown that carbon and oxygen in this gas are 
exceptionally overabundant compared to other elements, by some 400 times. We study the origin of 
this enrichment under two opposing hypothesis, preferential production, where the gas is produced with 
the observed unusual abundance (as may happen if gas is produced by photo-desorption from C / 0-rich 
icy grains) , and preferential depletion, where the gas evolves to the observed state from an original solar 
abundance (if outgassing occurs under high-speed collisions) under a number of dynamical processes. 
We include in our study the following processes: radiative blow-out of metallic elements, dynamical 
coupling between different species, and viscous accretion onto the star. We find that, if gas viscosity 
is sufficiently low (the conventional a parameter < 10"'^), differential blow-out dominates. While gas 
accumulates gradually in the disks, metallic elements subject to strong radiation forces, such as Na 
and Fe, deplete more quickly than C and O, naturally leading to the observed overabundance of C 
and O. On the other hand, if gas viscosity is high (a > 10~^, as expected for this largely ionized disk), 
gas is continuously produced and viscously accreted toward the star. This removal process does not 
discriminate between elements so the observed overabundance of C and O has to be explained by a 
preferential production that strongly favors C and O to other metallic elements. One such candidate 
is photo-desorption off the grains. We compare our calculation against all observed elements (^ 10) 
in the gas disk and find a mild preference for the second scenario, based on the abundance of Si 
alone. If true, (3 Pic should still be accreting at an observable rate, well after its primordial disk has 
disappeared. 

Subject headings: circumstellar matter planetary systems: formation planetary systems: protoplan- 
etary disks stars: individual (/? Pictoris) 



1. INTRODUCTION 



The young (~ 20Myr, 



Mentuch et"all I2008D . nearby 



e yc 

(19 pc, Ivan LeeuwenI [20071 ) main-sequence (A5 V) star 



/3 Pictoris has been known to harbour circumstellar 
(C S) gas eve r since its characterization as a shell star 
by ISlettebakI ([T975h With the discovery of CS dust 
by IRAS (|A umanni Jl985l) located in an edge-on disk 
(|Smith fc Terrilelll984D . the interest in the CS gas was 
renewed with more detailed observations of the circum- 
stellar absorption lines in both the optical and the ul- 
traviolet. The low estimated lifetime of CS dust im- 
plies that it is probably replenished on short timescales 
thro ugh destructive collisions, c haracteristic for a debris 
disk ([Backman fc Parescelll993[) . 

Early attempts to understand the properties and origin 
of the CS gas were hampered by the lack of constraints on 
what was believed to be the dominant gas components: 
H, C, N, and O. Observations of these elements are dif- 
ficult because transition lines are mostly in the far-UV, 
where the terrestrial atmosphere blocks radiation and the 
stellar photospheric emission is very weak, or in the in- 
frared (IR) to far-IR for thermal emission from molecular 
(H2) or hyper-fine transitions (O I, C II), where again 
spectroscopic observations from space are required. An- 
other problem was the unknown spatial distribution - 
absorption spectroscopy does not reveal where along the 



line of si ght the gas is located, making solutions degener- 
ate (e.g. iLagrange et al.lll995l ). In particular, the stellar 
radiation force on e.g. Na I and Ca II is known to exceed 
gravity by more than an order of magnitude, yet from 
the absorption line profile the gas is observed to be at 
rest with respect to the star . The problem was revisited 
when lOlofsson et al.l (|2001h spatially resolved emission 
from Na I in orbit around the star, showing the gas to 
be widespread in the disk. Something must prevent the 
gas from escaping on free-as cend trajectories, but limits 
on th e column density of H2 (iLecavelier des Etangs et "aTl 
I2001h together with limits on H I (jFreudling et al.lll995f ) 
shows that H cannot be dense enough i n the disk to sig- 
nifica ntly brake these metal species ijBrandeker et al.l 
[200I . To solve this enigma, a variety of bra king scenar- 
ios were studied by [Fernandez et aLl (120061 ). who found 
that an overabundance of C by a factor of at least 10 over 
solar abundance would be sufficient for charged ions in 
the disk to behave as a single fluid of low effective radia- 
tion pressure and remain in orbit. C is an effective brak- 
ing agent because of its high abundance, low absorption 
coefflcient in the optical (resulting in a negligible radia- 
tion force) , and an appreciable ionization fraction. Since 
the observed elements experiencing a high radiation force 
spend a tiny fraction of their time in the neutral state, 
it is sufflcicnt that only ions arc efficiently braked. Ob- 
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servations by the Far-UV Space Explorer (FUSE) sub- 
sequently showed C to indeed be overabundant by 20 x 
with respect to other species (iRoberge et al.|[2006D . Re- 
cently, detailed absorption spectroscopy by HARPS at 
the ESO 3.6 m telescope ([Brandekeri 120111) and far-IR 
spectroscopy by Herschel/PACS of C II 157.7 emis- 
sion from /3Pic (jBrandeker et al.ll2012[ ) have inferred an 
even higher overabundance of both C and O, up to 400 x 
the solar ratio. 

The large amounts of C in the CS gas explains why the 
metallic species are not removed by the radiation pres- 
sure, but raises the new que stion where this huge 4 00 x 
overabundance comes from. iFernandez et al.l ()2006f ) ar- 
gued that the CS gas cannot be primordial, that is, a 
remnant from the initial star-forming nebula. Instead, 
the gas is thought to be secondary, and is produced 
by the dust grains in the debris disk, either through 
evaporat ion during mutual colhsions (hereafter the ECG 
scenario. iCzechowski &: MannI [2007|) . or through photo- 
desorption by UV radiation (hereafter the PDG scenario, 
IChen et aLll2007[) .^ The two proposals produce gas with 
different chemical patterns: the former should lead to 
gas with solar-composition (hydrogen poor), the latter 
to C/0 rich gas. 

If the observed abundances reflect those at production, 
these observations directly reveal the production mech- 
anism, and more importantly, the chemical composition 
of the dust grains. The latter property is currently un- 
constrained: we do not yet know whether the grains are 
mostly icy or rocky, or if they are depleted in carbon like 
bodies in our Asteroid belt. In the /3Pic gas disk, not 
just C and O, but elements like Na, Mg, Al, Si, S, Ca, 
Cr, Mn, Fe, Ni are also observed. This affords us an un- 
precedented peek into the chemistry of planet formation 
in outer planetary systems. 

However, the current abundances of elements in the 
gas might not directly reflect those at production, as the 
abundances could have evolved in time. In particular, 
the various degrees of radiation force experienced by dif- 
ferent species may result in some elements being removed 
at higher rates than others. Since the radiation force on 
both C and O is negligible, while it is significant on most 
other observed elements, this could potentially explain 
why C and O appear enriched with respect to other ele- 
ments. In this paper, we investigate the differential de- 
pletion of various elements in the /3 Pic disk, in an effort 
to explain the enrichment of C and O and to constrain 
the chemical composition of the dust grains. In detail, 
our paper is organized as follows. We describe our model 
in S}2] with some of the detailed derivations delegated to 
Appendices A and B. The results are presented in H2.41 
Wc discuss in i^lthe implications of these results in the 
broad framework of planet formation and summarize our 
work in ^ 

2. MODEL 

We study the evolution of CS gas using a simple 
one-zone model. We consider elements from Li to Ni 
with the exceptions of noble gases. This gas is most 
likely second-generation gas, reproduced from solid ma- 

^ The evaporat ion of comet-like bodies falling into the star (e.g. 
IBeust et al.lll989l ) has also been suggested to be responsible for 
producing gas very close to the star. 



terials after the primordial disk has dispersed. As such, 
the element hydrogen should be very under abundant, 
together with helium and other noble gases that do 
not easily condense. We include only the neutral and 
first ionized states of the elements. Higher ionization 
states are negligible in the /3 Pic disk because ionizations 
are mostly caused by the stellar photospheric spectrum , 
which has few very har d photons ([Fernandez et al.ll2006l : 
IZagorovskv et aLll2010D . The element C, being abun- 
dant and significantly ionized, and experiencing little ra- 
diation pressure, is both the main donor of free electrons 
and the main braking agent (C II ions) for metal species 
that are radiatively accelerated. In contrast, the other 
abundant element, O, remains largely neutral and is con- 
sidered inert. We consider four processes: 1) Gas pro- 
duction, which we assume to be continuous and due to 
an unspecified mechanism. 2) Outward radiation drift. 
Different gas species experience the radiation force to 
different degrees, parametrized by the radiation force co- 
efficient /3 that is the ratio between the radiation and 
gravity force. For /3 > 1, radiation dominates over grav- 
ity. But because of the centrifugal force, /3 > 0.5 is suffi- 
cient for particles born on circular orbits to be unbound 
from the star. 3) Field particle interaction. Particles 
that are accelerated outward will interact with the field 
particles, C II, through Coulomb scattering. 4) Viscous 
accretion. Since the gas disk is significantly ionized, it is 
reasonable to assu me that it is subject to the magnetoro- 
tational instabilitv (|Balbus &: Hawlevll99ltlHawlev et al.l 
I1995D which allows the gas to viscously accrete into the 
star. In the following, we will investigate how these four 
processes compete with each other to determine the evo- 
lution of elemental abundances. 

We ignore the following processes: scattering between 
ions and neutral gas, and between ions and charged dust 
grains. Due to Coulomb focussing, ions like C II are 
several orders of magnitude more efficient as a braking 
agent, when compared to neutr als.'^ Dust, while consid- 
ered by IFernandez et al.l (|2006D as an effective braking 
agent, is much less efficient than C II when the number 
density of C II is as high as observed. 

The mass of the star is ^ 1.75 and the 
age of the system <ago ~ 20Myr. The radia- 
tion sp ectrum of /3 Pic is adopted from iHauschildt et al.l 
()1999() . and the radiatio n pressure is ca lc ulate d fol- 
lowing the pro cedure of IFernandez et al.l (|2006[ ) and 
IZagorovskv et al. (2010. ). All elements are assumed to be 
in pho to-ionization / recombin ation equilibrium. Accord- 
ing to iBrandeker et al.l (120041 ) , gas (at least the element 
sodium) follows a spatial distribution similar to that of 
the dust, with most gas mass located at r ~ 100 AU and 
extending outward to hundreds of AUs. We simplify 
this into a one-zone model with r ~ 100 AU, assuming: 
a disk temperature of T^isk ^ 50K (jZagorovskv et aLl 
120101 ) and a carbon density of iVc = 50 cm""^ , as is ap - 
propriate for the mid-plane at 100 AU ()Brandekeill20 1 It ) . 

2.1. Radiative Acceleration and Collisional Braking 

A particle (ion or neutral atom) with /3 > 0.5 is accel- 
erated outwards by the radiation force , accompanied by 

This applies even to braking atoms, as polarizability makes the 
ion-neutral cross-se ction several times larg er than that of neutral- 
neutral interaction (IFernandez et ani2006l) . 



Composition of the (3 Pic gas disk 



3 



□ 



neutral drift velocities 



A 

Mg 



AAA 
O AA^A 



Sep, a 



p_|^ ™_ _ _ Mnfe_ 



"Trianqte:" V 
Square: v^.^^j 
Cross: v„„,, 



10 15 20 

otomic number 



10* 



10° 



effective drift velocities 
P 



Si 

A A 
O 



Ti pp MrFe™ 



% A 



S AA 



□ 



Tnangle: Vj;j„*(T— 
Square: v„j.^j*f^ 
Cross: u. 



10 15 20 

otomic number 



Fig. 1. — Left panel: neutral drift velocities, Dncu,i (triangles), iincu,2 (squares) and Dncu (crosses), plotted against atomic numbers 
for elements with /3 > 0.5. Most neutral species satisfy ?;neu,2 < '^neu.ii or, drift velocities limited by the finite ionization time. Right 
panel: effective drift velocities nx (red crosses), and respective contributions from the neutral state (/xfncu) and from the first ionized 
state ((1 — /x)i'ion), for elements that satisfy /3 > 0.5 in one or both of these states. Drift in the neutral state dominate the effective drift, 
with the exception of elements Si and Ca. Results are based on A^c = 50 cm""^, r = 100 AU, T^isk = 50 K, M* = 1.75 Mq, and the 
radiation spectrum of a /3 Pic-like star. 



random scattering via collisions with field particles C II. 
Multiplying the radiative acceleration by the mean-free- 
time after which the particle suffers a strong (90 deg de- 
flection) scattering, we obtain its equilibrium outward 
drift velocity. Barring other processes, this velocity de- 
termines the gas removal rate from our one-zone model. 
For an ion "x" , this velocity is (see Appendix B) 
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where mc and mx are the masses of a C atom and the ion 
"x", respectively, r is the radial distance to the central 
star of mass M^,, A^cii is the number density of C II, and 
Tdisk is the temperature of the gas disk. 

For a neutral atom "x" , the equilibrium drift velocity 
is (see Appendix A) 



/ Tfi 
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where mp is the mass of the proton, and Px the polar- 
izability of neu t ral atom "x". We ad opt values from 
iJohnsonI ()2011h (http://cccbdb.nist.gov) and where un- 
available we simply take it to be Px = 10~^®m~3. 

The typical ionization time is short for neutrals in the 
/? Pic disk. This warrants considering the limiting ve- 
locity an atom can be accelerated to before it becomes 
ionized: 



GM^ 1 



^^ncu,2-^/3^f^- 1000/3 
1 



Fx 



AU 



cms ^, (3) 



where G is the gravity constant, Fx^au is the ionization 
rate of the neutr al atom "x" at r = 1 AU. Following 
iBrandekerl ()201lD 3. we adopt the following to be the 

^ Note that we use a different naming convention compared to 
IBrandekerl II2011I ). The two neutral limiting velocities Vneu.i and 



actual drift velocity of the atom, 

_ 7 ^nou,2 

7 + 1 



(4) 



where 7 = t;neu,l/?^neu,2- 

Since any given element switches between the neutral 
and the ionized state s, we combine Wneu and ''^ion to 
obtain the effective drift velocity for element "x" as 



= /xWnou + (1 - /x)wi( 



(5) 



where /x is the neutral fraction of element "x" that quan- 
tifies the fraction of time that is spent in the neutral 
state. 



ncu,li ^ncu,2; ^ncu 



and fin 



Fig. [T] shows the calculated d, 
for all elements with /3 > 0.5. We observe that: 



1. most elements' neutral drift velocities Wncu ai'G close 
to Unou,2 (the ionization limit), except for Be and 
P, which are less rapidly ionized. 

2. neutral drift (vncu/x) dominates the final effective 
outward drift for most elements, except for Si and 
Ca. The ion drift is almost completely suppressed 
by the strong coupling with C II. 

2.2. Radiative Acceleration versus Viscous Accretion 

The outward drift of radiatively accelerated metals is 
moderated by the inward diffusion in the viscous accre- 
tion disk. The latter arises in our partially ionized disk 
since the strong coupling between the magnetic field and 
the gas allows the onset of magnetorotational instability 
(jBalbus fc Hawle"vlll991[ ). This instability leads to angu- 
lar momentum transport t o which one might ascribe a 
a-like turbulent diffusivity (jShakura fc Sunvaevlll973[ ) 



(6) 



where Vg is the local sound speed, h the local scale height 
and r^K the Keplerian angular frequency. We set h = 
0.033AU(r/AU)^/* and estimate the typical timescale s 



'fncu,2 correspond to I'drift and Dion in lBrandekerl 1)20111 ). while vion 
in this paper denotes the limiting velocity of the ionized tracer 
particles. 
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for viscous accretion and radiation driven drift to be 
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100 cms 



100 AU 
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y^8) 



Note that ivis is the same for all elements, while ^drift is 
specific to element x. Equating ivis with tdrift we find that 
the critical a, at which the inward flow driven by viscous 
accretion is comparable to the outward flow driven by 
the radiation force, is given by 

acr ^ 0.3 (—^^^—^ . (9) 
V 10-^ cms ^/ 

As is shown in the right panel of Fig. [TJ many elements 
(Na included) have u^^ lO'^cms"^. So unless a <^ 
10^^, the inward viscous diffusion can compete against 
the outward radiative drift which differentially depletes 
element x. 

2.3. Evolution of the Abundances 

The key quantity out of this study is the abundance ra- 
tio. We define an abundance ratio between carbon and 
element x as one that is normalized to the solar value 
(|Anders fc Grevessd [Ml), [C/x] = logio(iVc/^x) - 
logio(iVc/A^x) Solar- The observed value is A^c/-^Na ~ 
400 (iVc/A^Na) Solar Or [C/Na] w 2.6 . It is the goal of this 
paper to understand how this large factor comes about. 

We study a model where the gas is continuously pro- 
duced at a constant rate, Sc, Sx, and with a fixed abun- 
dance ratio [C/x]src- The value of [C/x] in the CS gas 
evolve in time as the element x is subject to differential 
removal by radiation pressure while all elements viscously 
diffuse to the star at the same rate. The evolution is de- 
scribed by the following equations 



dN^ Nx 

dt ^vis 

which have solutions 

iVc = S'civis(l - e 



Nx 

idrift ' 



(10) 



iVx = ^xteff(l-e-*/*"«), (11) 

where <eff = ivis Adrift /(ivis + Adrift). Therefore, the evo- 
lution of the abundance ratio is 



[C/x]-[C/x],,,+logi 



cff 



1 



-*/tvis 



1 



-t/tctf 



(12) 



Elements that do not experience significant radiation 
pressure are easily incorporated into the above analysis 
by setting idrift —J- oo. 



2.4. 

Na, 



Results 



We first focus on Na, the element wit h most de- 
tailed observations (Brandeker et"aLl 120041 : iBrandekeil 
120111 : iBrandeker et al.ll20"T^ .^ Our calculations indicate 

* While Na in the /3 Pic disk mostly exists in ioni zed form that re- 
mains invisible to us, extensive modelling efforts IIFernandez et aLl 




Fig. 2. — Mapping the present abundance ratio [C/Na] (nor- 
malized to the solar abundance with logarithm scale) in the 
a — [C/Na]src plane. The contour line marked v^ith "obs" denotes 
the abundance ratio of carbon and sodium ([C/Na] = 2.6) that is 
about 400 times enhanced with respect to the solar abundance, as 
suggested by Brandeker ( 201ll). As expected, the present [C/Na] 
value increases with [C/Na]si.c and decreases for larger a. At small 
viscosity, e.g., a < 10~^, the viscous timescale is longer than the 
current age and [C/Na] becomes independent of a. In order 
to reproduce the observed [C/Na] value, the abundance of car- 
bon in the gas source should be either near solar with low viscosity 
{a < 10"'^) or super-solar with high viscosity (a > 10~^); a sub- 
solar abundance case can be excluded. Two crosses mark a typical 
solar (case I) and a typical super-solar cases (case II) that both 
reproduce the observed carbon enhancement. In Fig. 4, we study 
these cases in more detail. 



that there are two possibile scenario to reproduce the ob- 
served value of [C/Na] w 2.6. One is that the metallic gas 
is produced at the solar-composition, [C/Na]src ^ 0, but 
experiences a negligible turbulent viscosity, a < 10"'^. 
We call this 'preferential depletion . The other is that 
the gas is produced with enriched C (and O) abundances, 
[C/Na]src > 0, and it is advected into the star with larger 
a. We call this 'preferential production'. Our calculation 
places a constraint on the relationship between a and 
[C/Na]src, as shown in Fig. [2] But short of other the- 
oretical or observational evidences, both scenarios can 
explain the observed C/Na overabundance. 

Now we turn to examine other metallic species in the 
/3 Pic disk, in the hope of extracting further constraints 
beyond what is achieved using Na alone. We study 
two specific sets of parameters, both reproducing [C/Na] 
« 2.6 at current time: solar abundance production and 
a = 10~^ (preferential depletion); and super-solar pro- 
duction, [C/x]src — 2 and a = 0.1 (preferential produc- 
tion). We assume for both cases that [C/0] = 0, or Oxy- 



1200 6: Zaeor ovskv et aT]|2010l ) have led to estimate of the electron 
density (and therefore Na Il/Na I ratio) that consistently repro- 
duce the observed ionization fractions in various elements. Tliese 
are listed in Table [Tl 
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gen is produced in tandem with carbon. This is based on 
the considerations that O is similar to C in dynamics (ex- 
periencing only viscous diffusion), and that the observed 
[C/0] K. (jBrandeker et al.ll2012l ) . The resuUant [C/x] 
are plotted in FiglSl and summarized in Table [TJ for all 
elements between Li and Ni, with the exceptions of noble 
gases. 

In the case of solar production (top panel of Fig. [3]), 
the present day elemental abundances are determined 
by their effective drift velocities (right panel of Fig. [T]), 
with elements that experience strong radiation drift to be 
severely depleted relative to C and O. This is the case for 
Na, Mg,P, S, Cr, Mn, Fe, Ni.'^ But we expect elements 
that are not radiatively accelerated, like Al, and to some 
degree. Si, to be undepleted relative to C and O. This 
expectation contrasts strongly against the observations 
of Si, which is detected at a value 100 times lower than 
our model prediction.^ This difference is not explain- 
able by hiding most Si in higher ionization states: in the 
/? Pic disk. Si should be mostly in the first ionized state, 
which is eas ily observable. The ob servation of Al is more 
ambiguous. iLagrange et al.l (|1998[ ) only reported a lower 
limit to Al since the Al II line at 1670.79 A is very satu- 
rated. Based on the observed line profile, a rough back- 
of-the-envelope estimate yields that Al is likely ^ 10^ 
more abundant than the lower limit. If true, this would 
place the observed Al abundance in between values ex- 
pected for the two opposing scenario. The abundance 
of Al being crucial for deciphering the gas origin, a re- 
analysis of its line profile seems warranted. 

The agreement between theory and observation is 
mildly improved in the case where the gas production is 
super-solar (in C and O), and where gas removal is dom- 
inated by turbulent diffusion (bottom panel of Fig. [3]). 
The discrepancy, present still in Si, is a factor of 10 or 
smaller. For instance, for a production abundance of 
[C/x] = 3, the predicted Si is ^ 5 times higher than the 
observed value. In the following, we discuss the implica- 
tions of these results. 

3. DISCUSSION 

Here, we consider other astrophysical constraints on 
our scenario, and what future observations may bear on 
our current deliberation. The implications on the pro- 
cess of gas generation, as well as how our results are 
projected in the broad picture of planet formation, are 
also discussed here. 

3.1. Viscosity in the gas disk 

Since the gas disk has a high ionization fraction (^ 
20%), one expects it to be subject to the magnetorota- 
tional instability (MRI), which has been shown numer- 
ically to l ead to turbulent vi scosity of order a ^ 10^^ 
or higher (jHawlev et al.|[l995l ). We consider briefiy the 
complication due to the presence of charged dust grains 
in the debris disk. These grains couple s trongly to the 
charged particles ([Fernandez et al.l 120061 ) . The mean- 
free-time for electron-dust scattering is longer than the 
typical electron gyration time, so the electrons are ex- 
pected to be tightly coupled to the magnetic field. The 

^ The CS gas of /? Pic is likely of secondary origin, so a lower 
than solar abundance in H is not surprising. 

® This conclusion is not significantly altered when varying a. 



grains, however, may substantially enhan ce the resistiv - 
ity The latter effect is known to reduce a ( Balbusir2009D . 
down to values of orde r lO"** or lower (jFleming et al.l 
l2000tlBai fc Stonell201ll ). A more detailed study on MRI 
in the presence of charged dust is required to determine 
the realistic a value. 

3.2. The rate of gas production and accretion 

We obtain the current total mass of the gas disk by 
combining the observed column densities of C in Table 
[T] with an assumed scale height h ^ O.lr. Further as- 
suming that C and O dominate the total gas mass and 
are produced at the solar ratio relative to each other, we 
estimate the required gas production rate to be 



I 1.2x 10i3(^) gs-Mfa>7x 10-4 



dt 



X 10^" gs 



10 „c.-i 



(13) 



otherwise. 



We compare these values against various models for 
gas production. For the solar abundance case, 

iCzechowski &: MannI (120071 ) estimated the gas produc- 
tion rate in their ECG model as 5 x 10^^ gs^^ for silicate 
grains and 2 x 10^"^ gs"^ for icy grains, in the /3Pic disk. 
These are roughly consistent with our above requirement. 

Photo-desorption is a class of model where molecules 
and atoms are e jected from grains by UV photons. 
IChen et all (|2007D estimated the photo-desorption rate 
for Na, finding a Na production rate of ~ 1.33 x 10^"^ 
for /3Pic. This corresponds to a total gas production 
rate of 1.0 x 10^^ gs"^ if photo-desorption efficiency is 
the same for C, N and O, as for Na. lOberg et all (|2007[ ) 
performed experiments of photo-desorption of CO and 
CO2 gas from icy grains. They found rates of order 10"^ 
molecules per UV photon, or, 100 times higher than, 
e.g., those used by IChen et"aLl (|2007D for Na. So ff the 
dust grains have solar abundance and if all C and O are 
stored in the form of CO or CO2 ice, photo-desorption 
would naturally yield a gas that is enriched in C and O by 
~ 100 times relative to Na (and perhaps other elements). 
Moreover, the gas production rate would be ~ 10^^ g s"^. 
So interestingly, both classes of models have no difficulty 
in delivering the required amount of gas. We cannot ex- 
clude either model using this discriminant. 

The dust disk is also losing grains through grinding 
down and radiation pressure push-out. Taking a vertical 
optical depth of 10^'^ for the grains at 5 /im size (blow- 
out size), we arrive at a ra te of ^ lO^'^gs"^ (also see 
iThebault fc Augereaul [20071 ) . So for high a {a > 0.1), 
the efficiency of gas production is almost comparable to 
the dust removal rate. We then have the intriguing pos- 
sibility that dust in the (3 Pic disk may not be all blown 
away but is partially accreted back to the star in the 
form of gas. This possi bility was also considered by 
IZuckerman fc So"n3 (|2012[ ) for the bright debris disk of 
49 Ceti. 

In the super-solar case, we expect to see the gas be- 
ing predominantly channelled into the central star. Such 
an accretion could potentially provide an explanation for 
the X-ray emission observed from /? Pic: the required ac- 
cretion rate to ex plain the X-ray emi ssion is 6.6 x lO^'^ 
- 6.6 X lO^^gs-i (jHempel et al.ll2005D . compatible with 
the above production rates for high a's. This gas, more- 
over, should be highly enriched in C and O and relatively 
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case I: solar abundance, preferential depletion 
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cose II: super-solar abundance, preferential production 
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Fig. 3. — Abundance pattern in the /3Pic disk, normalized to the solar value, and with all elements scaled by the element Na. The 
blue boxes, arrows, and errorbars represent the observed values and the red crosses the results of our calculations with the vertical red 
errorbars indicating effects of one source of theoretical uncertainties, when we vary the ionization fraction by a factor of 3 each way from 
those calculated in table 1. The top panel shows the case when all elements are produced at the solar pattern, and the viscosity is low 
(a = 10~*); the bottom panel shows the case when all elements are produced at the solar pattern except for C and O which are prod uced 
100 times more abundantly, and when a = 0.1. The size of the blue box reflects a typical uncertainty of 40% l ILagrange et al.lll998l) . and 
the arrows indicate upper limits (for P) and lower limits (for Mg and Al which have saturated absorption lines), respectively. Calculated 
values for elements that are as yet undetected are plotted as small black crosses. Also see Table [T] for tabulated values, as well as references 
for the observational ones. 
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TABLE 1 

Column densities and Abundunces: observation vs. theory" 



atomic 


element 


neutral column 


ion column 


fx 






[x/Na] 


[x/Na] 


References 


number 




density 
[cm-2]6 


density 
[cm-^Jf- 


(obs)= 






(obs) 


(cal)'* 




6 


C 


1.4 X IQi' 


2.9 X 10^'' 


0.33 




0.52 


2.58 


2.54 


2.55 


[1] 


8 


O 


8.6 X 10" 








~ 1 


2.58 


2.54 


2.55 


[1] 


11 


Na 


(3.4 ±0.4) X lOi" 


[9.2 X 10^2] 




3 


7 X 10-3 











[2] 


12 


Mg 


2.5 X 10" 


> 2.4 X 10" 


< 1 X 10-2 


8 


6 X 10-3 


> -0.83 


0.15 


0.11 


[3], [4] 


13 


Al 


< 2.7 X 10^2 


> 4.4 X 10^2 


< 0.4 


1 


3 X lO-'^ 


> -0.49 


2.54 


0.55 


[3] 


14 


Si 


< 1.1 X 10^^ 


1.1 X 1014 


< 1 X 10-1 


1 


1 X 10-3 


-0.16 


1.62 


0.53 


[3] 


15 


P 


< 7 X 10" 


< 9.2 X 1013 






0.22 


< 1.8 


-0.83 


-0.71 


[5] 


16 


S 


5.4 X 10^2 


[4.9 X 10i3] 






0.1 


-0.12 


0.13 


0.09 


[3] 


20 


Ca 


< 2 X 10^ 


2.6t^;^ X 1013 


< 1.4 X 10-1 


5 


7 X 10-"' 


041+0.6 


1.18 


0.49 


[3], [6], [7] 


24 


Cr 


[3.5 X lO^"] 


4.8 X 10i2 




7 


4 X 10-3 


0.37 


0.18 


0.12 


[3] 


25 


Mn 


< 2.4 X 10^0 


3.8 X 1012 


< 6 X 10-3 


9 


4 X 10-3 


0.41 


0.47 


0.28 


[3] 


26 


Fe 


4.9 X 10" 


1.6 X 1014 


3.1 X 10-3 


9 


8 X 10-3 


0.44 


0.46 


0.28 


[3], [7] 


28 


Ni 


< 7.6 X 10^0 


1.5 X 10i3 


< 5 X 10-3 


5 


4 X 10-3 


0.29 


0.91 


0.43 


[3], [7] 


References: 


[11 Brandeker (2011) [21 Brandeker 


et al. (2004). [31 L, 


igrange et al. (1998), 


[41 Vidal-Madiar et al. 


1994) 








ISlRoberee 


et al. C2006'). [61 Vidal- 


Madiar et al. C1986), [71 Crawford et 


al. 


(1994) 











Note. — ° Only elements that have reported column densities (for either neutral or first ionized state) are presented here. 

Note. — Observed values for column densitie s from absorption line studies (except C/O/Na where emission lines are detected). Most 
error bars are of order 40% of the measured values I ILagrange et al.lll998l ). unless explicitly listed. Where multiple absorption features from 
the same lower state are present for the same species, we adopt the largest value of column density. If the lower states are different, we 
add up the column densities. Values in square brackets are calculated assuming ionization equilibrium {fx cal). 

Note. — Neutral fraction 

Note. — Abundances relative to solar, normalized by Na. Values in front of the semicolons are for the solar case, while those behind 
are for the super solar case (top and bottom panels of Fig. [S] respectively). 
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depleted in other metals. This prediction can be tested 
using data like those from lBouret et all (|2002[) . 

3.3. Carbon and Oxygen Ice-Lines 

Our result that C and O have to be produced at solar 
or even super-solar abundance with respect to other ele- 
ments has interesting implications for the environments 
of proto-planetary disks. If the gas originates from grains 
~ 100 AU away from /3Pic, then C and O would have to 
be of at least solar abundance in these grains. This con- 
trasts with both the asteroid belt bodies and the more 
distant cometary bodies in the solar system. 

In the proto-solar nebula, the ice-line for water ice (of 
condensation temperature 135 K) is calculated to lie at 
~ 1 AU, while that for CO ice (of condensation tem- 
perature 20 K) is located at 10-100 AU, depending on 
disk geometry (jChiang fc Goldreichlll997| ). The flatter 
the disk geometry, the closer in the ice line. Solid bod- 
ies formed between the two ice-lines would have been 
rich in water but carbon poor, with densities close to 
Igcm^'^. This is indeed observed for the satehites o f 
Jupiter and Saturn (see review bv IWong et all l2008l) 
and for smaller Kuiper belt objects (jGrundv et al.ll2007l : 
iStansberrv et al.l 120061 ).'' Comets are also known to be 
significantly depleted in carbon relative to water ice, in- 
dependent of whet her they arrive from the K uiper belt 
or the Oort cloud (iMumma fc Charnlevl[20lTh . This in- 
dicates that the carbon ice-line lies outside ~ 40 AU in 
the solar system. 

Scaling from the solar model, with a stellar luminosity 
of 8.7Lq for /3Pic, the ice-line for CO should have lied 
around 30-300 AU. Here, our study indicate that the car- 
bon ice-line should have lied inward of ^ 100 AU, or, the 
primordial disk is not very flared. 

3.4. Future Refinement 

In the present paper, we have described a local model 
which assumes that gas is produced at 100 AU, where 
most of the dust is located, and compute how the ra- 
diation force and accretion differentially act on the gas 
to evolve the abundances. In reality, gas is observed 
to extend from at least 13 AU from the central star 
([Brandcker et al. 2004 ) to hundreds of AU. If the gas 
is produced where most of the dust is located, the ques- 
tion becomes how the gas is transported inwards. One 
possibility is accretion triggered by MRI. To model this, 
however, requires a stratified model. Another assump- 
tion of our analytical model is that C II braking is ef- 
ficient already from start, resulting in a constant drift 
velocity and a constant production rate for all el- 
ements "x" . A more realistic scenario would be that the 
gas builds up over time with C II braking becoming ef- 
ficient only after the C density has reached a sufficient 
level. A more refined model should thus take into ac- 
count both the spatial extension and the evolution of 
the gas, to provide more precise abundance predictions 
for various elements, yielding stronger constraints on the 
origin of the gas around (3 Pic. 

4. SUMMARY 

Surprisingly, large Kuiper-belt objects like Pluto and Eris 
are too dense (~2g cm~^) to be water enriched IINuU et al.lll993l : 
ISicardv et aLlljoTiri . 



In this paper, we investigate the unusual elemental 
abundance of the /? Pic gas disk. We find that to attain 
the observed C/0 rich gas using our model, the C and O 
gas must be produced at either solar or super-solar abun- 
dance compared to other elements; the sub-solar case is 
excluded. In the case of solar abundance, a low gas vis- 
cosity of a < 10~^ is required, as well as a gas produc- 
tion rate of ~ l O^^gs"^. Such a rate is achievable in 
the ECG model (jCzechowski fc Mannll2007[) where gas is 
produced by high velocity micro-meteorites bombarding 
and sublimating the debris disk grains. In this scenario, 
the currently observed overabundance of C and O results 
from preferential depletion of other metallic elements by 
radiative acceleration. 

In the case of super-solar production of C and O, a 
high gas viscosity is required (a ~ 0.1), as well as a high 
production rate of lO^'^gs"^. Such a rate, as well as 
the exotic production pattern, can be satisfied by photo- 
desorption if the recent study on the photo-desorption of 
CO stands (jOberg et al.ll2007[ ). The mass in the /3Pic 
dust disk may be partially ground down to micro mete- 
orites that are pushed out by radiation pressure, partly 
outgassed and accreted into the star. The latter part may 
be significant. As this gas accretes into the central star, 
it may provide a natural explanation for the observed 
x-ray emission from (3 Pic, unusual for an A-star. 

Detailed comparison between our model and the ob- 
served abundance pattern leads us to prefer the super- 
solar case, but only by a small margin. Our preference 
is informed by the fact that the column density of Si is 
closer to (albeit still 10 times larger than) the observed 
value in the super-solar case. Moreover, we argue that 
the gas disk should support the magnctorotational insta- 
bility which yields a healthy viscosity. We propose the 
following tests to determine the dominant gas production 
mechanism: 

• A re-analysis of the Al II absorption data of 

grange Since Al should not be pref- 

erentially depleted by radiation pressure, its abun- 
dance ratio relative to C and O will inform us of 
the production pattern. 

• Analysis of UV a nd x -ray data like those in 
iBouret et all ([200l and iHempel et all (|2005l) . If 

the x-ray does arise from gas accretion, as opposed 
to chromospheric activity, we can use it to deter- 
mine the gas accretion rate, as well as the chemical 
abundances. 

• The radial profile of the gas disk. If the outgassing 
occurs at ^ 100 AU and most of the gas accretes 
into the star, we expect to see gas at smaller radii. 
This is a particularly exciting aspect beca use the 
giant planet, /3Pic b (jLagrange et al.ir2010( ). is be- 
lieved to orbit at ~ 10 AU and may interact with 
the gas in interesting ways. 

Otu' current model is limited in that the gas disk is 
approximated by a single zone at steady state. This, 
though unlikely to affect our basic conclusions, would 
need improvements to allow a more detailed comparison 
with observations. 
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APPENDIX 

A. EQUILIBRIUM VELOCITY DUE TO NEUTRAL-ION COLLISIONS: Vneu,i 

Following iBeust et al.l ()19890 . who studied the braking of ions by a neutral gas, we here consider a neutral atom 
braked by an ion gas (consisting of C II at a number density of iVcii). The average net momentum loss of the neutral 
atom in one collision is — mcv, where ttiq is the mass of one C atom and v is the neutral drift velocity relative to 
C II, with V its magnitude. On average, the collisions arc then equivalent to an effective force on the atom that can 
be expressed as 

V 



-NcuTrb^imcw — —k — v. (Al) 

Vcl 



Here 



/ 1 4e2p^ 



1/4 



V = 7 ^ (A2) 



is the largest impact parameter that can lead to a physical collision between the neutral and the ion, where eg is the 
permittivity of free space, e is the charge of an electron, is the polarizability of the neutral "x" , ^ ~ m^mc / {m^+mc) 
is the reduced mass, v^y is the neutral-ion collision velocity, and where 



4e2p, 7V2jj 

k = TTmc\ — (A3) 

V 47reo ^ 

Two parts are contributing to the coUisional velocity, 

VA-v^ + v, (A4) 
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where Vg is the sound speed of C II gas. The atom, on the other hand, is subject to the radiation force from the central 
star, which is 

Frad ^ P ■ (A5) 

Equating Fneu to Frad, the equihbrium velocity can be solved as 



_ fiad + v/f;:ld + 4fc^;.F,ad 

fnou,l — ■ 

This is the limiting drift velocity of a radiatively driven atom braked by neutral-ion collisions. If w S> (an assumption 
that here holds for most elements), then equation IA4I is reduced to v^i ~ w, leading to 



^ncu, 1 

TT \mcj \ J \ 47reo 



\ 1/2 / AT \ -1 , . 9 / n r \ / u \ -1/2 



mc J \mp{m^ + mc) ) V100cm-3y VlOOAU/ \Mq J 
which is the case considered by iBeust et al.l ()1989( ) . 



cms \ (A7) 



B. EQUILIBRIUM VELOCITY DUE TO ION-ION COLLISIONS: Vion 

Here we study the case of a tracer ion moving in C II gas. Following IBeust et al.l ()1989D . according to the classical 
theory of Coulomb scattering (here referred to as "collision"), the average net momentum loss of a tracer ion colliding 
with a C ion is 



(5p = — mcv(cos X — 1), (Bl) 



where x is the deflection angle given by 



47reo /i bvl 



tan^ = 7r— I-T' (B2) 



where h is the impact parameter, v^x is the ion-ion impact velocity as given by equation IA41 In IBeust et al.l (jl989() . 
they implicitly ignore Vs and take Vc\ v. In the case we consider here, however, the tracer ion would be efficiently 
braked by the C II, implying v <^ Vg and Vd ^ Vg. 
The equivalent force due to ion-ion collisions can be expressed as 



:uiVcii / {Sp)2Trb db 
Jo 



2Ncin:bl In ( -^-^ + 1 ) !!^mcvv, if 



^mcuTTbl In (^^^ (^1 + mcvv. (B3) 

Here 6max is the maximum acceptable value of the impact parameter, which can be taken as the Debye length 



Ad ~ \ ~ W £o|B7di^^ .j: ^ ^^^^ ^ T^is)^, (B4) 

V e^rie V e^A'cii 

where fee is the Boltzmann constant, Tg and Tdisk are the temperatures of the electrons and the disk, respectively, and 
Ue is the number density of electrons. If we assume C to be the major donator of electrons, then we expect Ue Nqii- 
bii is a characteristic impact parameter for the field ion- ion collision (obtained by solving equation IB2I for x = 90° and 
substituting ^ with toc): 

1 e2 1 1 e2 



thus 

^=4. .3/2 ,-34/2^,3/. ^^-1/2. 
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Equating Fio„ with F^ad, wc find the equilibrium velocity to be 



ff 



/3 



f 



41n(AD/6ii)/ \Trbf-rNciiJ \mc{mc + m^) 



1/2 



1/2 



'3.7/3 



mc{mc +TOx) 



1/2 



V fOOcm- 



(fOOAu) 



V fOOK 



3k \ -1 

cm s . 



(B7) 



